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Silicon  and  germanium  are  both  recognized  as  a  promising  anode  material  for  high-energy  lithium-ion 
batteries.  Si  is  best  known  for  its  superior  energy  storage  capacity,  while  Ge  exhibits  better  rate  capability 
and  cycleability.  To  better  understand  the  underlying  reasons  behind  their  lithiation  behavior  differences, 
particularly  the  enhanced  Li  transport  in  Ge,  we  examine  and  compare  Li-host  lattice  interactions  and 
dynamics  using  density  functional  theory  calculations.  At  the  onset  of  lithiation,  an  isolated  Li  interstitial 
is  found  to  form  polar  covalent  bonds  with  four  nearest  host  atoms,  while  the  degree  of  covalency  is 
noticeably  greater  for  Li— Si  than  Li— Ge  bonds.  The  relatively  stronger  Li— Si  interaction,  along  with  the 
stiffer  Si  lattice  tend  to  be  responsible  for  the  suppressed  Li  mobility  (Du  =  10-13  cm2  s-1)  in  c-Si,  as 
compared  to  the  c-Ge  case  (Du  =  10-11  cm2  s-1).  With  continued  lithiation,  Du  in  a-LixSi  increases 
significantly  from  10-12  to  1CT7  cm2  s-1  (x  =  0.14—3.57);  contrarily,  Du  in  a-LixGe  is  around  10-7  cm2  s-1 
and  less  concentration  dependent.  Our  analysis  shows  that  the  rapid  Li  diffusion  in  a-LixGe  is  directly 
related  to  the  facile  atomic  rearrangements  of  host  Ge  atoms  even  at  the  early  stages  of  lithiation. 
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1.  Introduction 
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Li-ion  batteries  (LIBs)  have  received  tremendous  attention  as 
they  power  a  wide  range  of  applications  from  portable  devices, 
electric  vehicles  to  various  renewable  energy  systems  [1—3]. 
Currently,  the  most  adopted  anode  material  is  graphite,  which  has 
good  cycleability  but  the  dendrite  formation  raises  safety  concerns 
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and  the  capacity  is  rather  limited  (372  mAh  g-1)  especially  at  high 
charge/discharge  rates.  Therefore,  in  order  to  satisfy  the  ever- 
increasing  energy  density/power  capability  requirements  and 
stringent  safety  standards,  there  is  an  imminent  need  to  find  new 
electrode  materials  with  superior  lithiation  properties.  Among  the 
alternatives  considered,  Si  stands  out  the  most  because  of  its 
impressive  capacity  (4200  mAh  g1  for  L^Sis  [4,5]),  safe  thermo¬ 
dynamic  potential  and  abundance.  Second  only  to  Si,  Ge  has  a 
relatively  high  theoretical  capacity  of  1624  mAh  g_1  (L^Ges  [6]),  a 
higher  electrical  conductivity  compared  to  Si  [7],  and  a  superior 
rate  capability,  up  to  1000  C  (full  lithiation  in  1/1000  h)  [6].  How¬ 
ever,  the  understanding  and  development  of  Ge-based  anodes  have 
gained  much  less  attention  likely  because  of  its  higher  price  relative 
to  Si. 

Being  in  the  same  column  in  the  periodic  table,  Si  and  Ge  share 
many  similarities,  including  the  disadvantages  of  undergoing  large 
structural  changes  and  volume  expansion  upon  lithiation,  which 
can  consequently  lead  to  early  capacity  fading.  To  overcome  this 
drawback,  many  ongoing  studies  have  focused  on  utilizing  Si  and 
Ge  of  different  forms,  such  as  thin  films  [8-12],  nanoparticles 
[13,14],  nanowires  [15-18],  and  alloys/composites  with  active/ 
inactive  elements  [19-22].  In  both  Si  and  Ge  cases,  nanostructuring 
seems  to  have  positive  impacts  on  enhancing  the  rate  capability 
and  reducing/preventing  electrode  pulverization,  thereby 
improving  the  cycleability.  In  comparison  to  Si,  Ge  of  comparable 
nano-architecture  is  able  to  withstand  much  faster  charging  rates 
with  noticeably  less  crack  formation  [23,24].  Furthermore,  there 
appears  to  be  subtle  differences  in  their  responses  to  electro¬ 
chemical  lithiation/delithiation,  as  demonstrated  by  recent  in-situ 
characterization  [24  .On  the  theoretical  side,  there  have  been  many 
studies  employing  density  functional  theory  (DFT)  to  examine  Li 
incorporation  in  Si  (crystalline/amorphous  bulks  [25-27]  and 
nanowires  [28,29])  and  a  few  on  Ge  [30-32].  Nonetheless,  the 
fundamental  understanding  regarding  the  nature  and  origin  of 
their  dissimilar  responses  to  lithiation  is  still  limited;  to  the  best  of 
our  knowledge,  no  atomistic  study  has  been  reported  to  investigate 
the  likely  overlooked  differences  between  Si  and  Ge  as  anode  ma¬ 
terial,  especially  regarding  their  lithiation  dynamics. 

In  this  paper,  on  the  basis  of  DFT  calculations,  we  examine  how 
Li  diffusion  is  affected  by  its  interaction  with  the  pure  Si  and  Ge 
matrices,  analyze  the  dynamic  behaviors  of  Li  as  well  as  the  host 
lattice  atoms,  and  look  into  the  impacts  of  Ge-alloying  on  anode 
performance.  The  fundamental  findings  explain  the  origin  of  the 
lithiation  behavior  differences  between  Si  and  Ge,  particularly  the 
significantly  enhanced  Li  transport  in  Ge,  and  thus  assist  the 
rational  design  of  the  next-generation  high  performance  Si-  and 
Ge-based  anodes. 

2.  Computational  methods 

The  calculations  reported  herein  were  performed  on  the  basis  of 
density  functional  theory  (DFT)  within  the  generalized  gradient 
approximation  (GGA-PW91)  [33],  as  implemented  in  the  Vienna 
Ab-initio  Simulation  Package  (VASP)  [34-36].  Spin  polarization  of 
the  Li-Si  (Ge)  system  was  also  examined,  but  appears  to  be  un¬ 
important.  The  projector  augmented  wave  (PAW)  method  with  a 
planewave  basis  set  was  employed  to  describe  the  interaction  be¬ 
tween  ion  cores  and  valence  electrons.  The  PAW  method  is,  in 
principle,  an  all-electron  frozen-core  approach  that  considers  exact 
valence  wave  functions.  Valence  configurations  employed  are  as 
follows:  ls22s1  for  Li,  3s23p2  for  Si  and  4s24p2  for  Ge.  An  energy 
cutoff  of  350  eV  was  applied  for  the  planewave  expansion  of  the 
electronic  eigenfunctions.  The  crystalline  Si  (Ge)  host  was  modeled 
using  a  216-atom  supercell  with  a  fixed  lattice  constant  of  5.457 
(5.777)  A;  the  effect  of  volume  relaxation  was  also  checked,  and 


turns  out  to  be  unimportant  as  the  216-atom  supercell  is  large 
enough  to  accommodate  one  Li  atom  with  no  significant  volume 
change  (less  than  1%).  For  geometry  optimization,  all  atoms  were 
fully  relaxed  using  the  conjugate  gradient  method  until  residual 
forces  on  constituent  atoms  become  smaller  than  5  x  10-2  eV  A-1, 
and  a  (2  x  2  x  2)  k-point  mesh  in  the  scheme  of  Monkhorst-Pack 
was  used  for  the  Brillouin  zone  sampling  [37].  Diffusion  pathways 
and  barriers  were  determined  using  the  climbing-image  nudged 
elastic  band  method  with  eight  intermediate  images  for  each 
hopping  step. 

The  model  structures  of  amorphous  a-LixM  alloys  (M  =  Ge  and 
Sii_yGey)  were  created  using  ab  initio  molecular  dynamics  (AIMD) 
simulations  based  on  the  atomic  configurations  of  a-LixSi  alloys  that 
were  previously  generated  using  the  combined  modified 
embedded  atom  method  (MEAM)  and  AIMD  simulations  (see  Refs. 
[31,38]  for  detailed  computational  methods).  The  interaction  be¬ 
tween  Li  and  Ge  is  very  similar  to  that  with  Si,  and  their  lithiated 
phases  tend  to  share  many  structural  similarities  [31,40,41]; 
therefore,  the  a-LixSi  structure  is  likely  a  good  starting  configura¬ 
tion  for  the  a-LixM  structure.  The  Si  atoms  (in  a-LixSi)  were  replaced 
by  Ge  atoms  accordingly  to  achieve  desired  composition  ratios;  the 
replacement  sites  were  carefully  chosen  to  ensure  homogeneous 
atomic  mixing.  The  model  structures,  each  containing  128  atoms, 
were  then  annealed  at  300  I<  for  1.5  ps  with  a  time  step  of  1  fs  to 
allow  sufficient  atomic  rearrangement  (the  annealing  temperature 
was  controlled  via  velocity  rescaling),  followed  by  geometry  opti¬ 
mization  using  a  (2  x  2  x  1)  k-point  mesh.  Periodic  boundary 
conditions  were  employed  in  all  three  directions,  and  for  each 
composition  three  independent  samples  were  considered.  This 
approach  can  provide  reasonable  Li-M  amorphous  structures  at 
significantly  reduced  computational  burden  compared  to  starting 
with  crystalline  initial  configurations.  Finally,  for  diffusivity  calcu¬ 
lations,  AIMD  simulations  of  8  ps  duration  with  a  time  step  of  1  fs 
were  carried  out  at  appropriate  temperatures  controlled  via  Nose- 
Hoover  thermostat. 

3.  Results  and  discussion 

We  first  examined  how  the  room-temperature  diffusivity  of  Li 
varies  with  Li  content  (x)  in  a-LixSi  alloys  using  AIMD  simulations 
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Fig.  1.  Predicted  diffusivity  of  Li  (DLi)  in  a-LixSi  and  a-LixGe  using  AIMD  simulations  at 
298  K.  The  predicted  DLi  values  for  single  Li  diffusion  in  c-Si  [42]  and  c-Ge  [31]  are 
indicated  by  red  and  blue  arrows,  respectively.  For  comparison,  DLi  values  calculated  by 
MD  simulations  with  an  embedded  atom  method  interatomic  potential  are  indicated 
by  *  [45].  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader 
is  referred  to  the  web  version  of  this  article.) 
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[Fig.  1].  Here,  a-LixSi  alloys  were  considered,  instead  of  their  crys¬ 
talline  counterparts,  because  Si  lithiated  beyond  the  first  charge 
cycle  is  most  likely  to  remain  in  the  amorphous  state  due  to  the 
sizable  kinetic  barrier  for  recrystallization  at  room  temperature 
[8,39].  For  each  alloy,  three  samples  are  averaged  to  calculate  the 
mean-square  displacements  (MSD)  of  Li  atoms  at  each  tempera¬ 
ture;  MSD  =  \Ri(t)  -  l?i(0)|2,  where  Rft)  is  the  position  of  atom  i  at 
time  t.  Based  on  the  MSD  profiles,  Dy  values  are  obtained  using  the 
Einstein  relation,  D  =  <MSD>/6t;  the  angular  bracket  denotes 
ensemble  average  over  the  AIMD  interval.  The  MD  duration  of  8  ps 
appears  to  be  sufficient  to  obtain  well-converged  results;  dis¬ 
regarding  the  first  2  ps,  linear  fits  over  a  time  interval  of  the 
following  6  ps  yield  the  Dy  values  [39].  Using  these  calculated  Dy 
values  at  different  temperatures,  an  Arrhenius  plot  of  ln(Dy)  versus 
1000 IT  is  constructed  based  on  D  =  Do  exp(-Ea//<T)  as  shown  in 
Fig.  SI ;  the  estimated  prefactor  D0,  diffusion  barrier  Ea,  and  room- 
temperature  diffusivity  Dy  values  are  summarized  in  Table  SI. 

Several  observations  can  be  made  based  on  the  above  calculation 
results.  Firstly,  in  comparison  to  Li  diffusion  in  c-Si  (where  Ea  and  Do 
are  predicted  to  be  0.62  eV  and  3.18  x  10'3  cm2  s'1,  respectively 
[42]),  the  Ea  in  a-LixSi  is  relatively  smaller  and  becomes  progressively 
more  so  with  increasing  Li  contents.  The  reduction  of  Ea  indicates  Li 
can  diffuse  more  easily  in  a-LixSi  of  larger  x,  which  is  not  surprising 
given  that  the  Si  host  matrix  would  be  softened  and  undergo  sig¬ 
nificant  disintegration  during  lithiation  [39,43].  Secondly,  contrary 
to  highly  concentration-sensitive  Ea,  D0  tends  to  be  less  affected  by 
compositional  changes  and  remain  relatively  constant  around 
10'3  cm2  s'1.  Note  that  our  predicted  D0  values  are  indeed  compa¬ 
rable  with  the  prediction  based  on  harmonic  transition  state  theory; 
according  to  which,  Do  =  vo a2  exp(-AS//<),  assuming  AS  (the  entropy 
difference  between  the  diffusing  atom  at  the  saddle  point  and  the 
initial  equilibrium  state)  is  close  to  zero,  vq  (the  attempt 
frequency)  ~  1013  s'1  and  a  (the  distance  between  adjacent  hopping 
sites)  ~  1CT8  cm,  D0  would  be  on  the  order  of  10'3  cm2  s'1  [44]. 
Lastly,  manifested  by  the  virtually  unchanging  D0  and  the  decreasing 
Ea  in  the  exponential  term,  Dy  in  a-LixSi  can  rise  by  orders  of 
magnitude  with  increasing  x,  from  xlO'12  cm2  s'1  (x  =  0.14) 
to  xlO'7  cm2  s'1  (x  =  3.57).  Such  trend  has  also  been  demonstrated 
by  recent  MD  simulations  although  some  discrepancies  may  appear 
due  to  different  choices  of  interatomic  potentials  [45]. 

Next,  same  analyses  are  carried  out  to  examine  Li  diffusivity  in  a- 
LixGe  alloys  of  varying  Li  contents  [Fig.  1  ].  In  c-Ge,  taking  Ea  =  0.44  eV 
(as  predicted  from  our  previous  work  [31])  and  assuming 
D0  =  10'3  cm2  s'1,  the  estimated  Dy  (-10'11  cm2  s'1)  at  room 
temperature  is  about  102  times  larger  than  that  in  c-Si.  Unlike  the  Si 
case,  Dy  in  Ge  is  found  to  increase  and  plateau  rapidly  with  lith¬ 
iation;  interestingly,  the  predicted  DLi  values  in  a-Li0.i4Ge  and  a- 
Li3.57Ge  are  of  the  same  order  (x  10  7  cm2  s'1).  Our  results  suggest 
that,  in  the  lithiated  Ge  matrices,  Dy  tends  to  be  less  concentration- 
sensitive  (as  compared  to  the  Si  case)  while  Li  diffusion  can  be 
extremely  facile  even  during  the  early  stages  of  lithiation.  This  dif¬ 
ference  in  lithiation  behavior  between  a-LixGe  and  a-LixSi  becomes 
progressively  smaller  with  increasing  Li  contents,  and  the  predicted 
Dy  values  are  both  on  the  order  of  10'7  cm2  s'1  when  x  =  3.57.  In  the 
following  sessions,  to  better  understand  the  origin  of  the  significant 
enhancement  of  Li  diffusivity  in  Ge,  particularly  at  the  early  stages  of 
lithiation,  we  analyzed  the  interaction  of  Li  with  the  host  matrix  as 
well  as  the  dynamic  behavior  of  host  lattice  atoms. 

3.1.  The  interaction  of  Li  with  c-Si  and  c-Ge  host  matrices  at  0  I< 

In  both  c-Si  and  c-Ge,  the  tetrahedral  (T)  site  has  been  identified 
to  be  energetically  the  most  favorable  interstitial  site  for  a  Li  atom, 
which  may  undergo  migration  by  jumping  between  adjacent  T- 
sites  via  the  hexagonal  (H)  transition  site  [42].  Despite  the 


structural  similarities,  the  predicted  diffusion  barrier  (Ea  =  0.44  eV) 
in  c-Ge  is  considerably  lower  than  0.62  eV  in  c-Si.  The  difference  in 
Ea  could  be  attributed  to  two  factors;  (i)  the  flexibility  of  the  host 
lattice  to  expand  and  allow  Li  passing  though  [31],  and  (ii)  the 
interaction  between  the  Li  and  host  atoms,  i.e.,  a  greater  Ea  would 
be  expected  if  the  diffusing  Li  interstitial  interacts  more  strongly 
with  neighboring  host  atoms  in  the  minimum-energy  T  state 
(relative  to  the  transition  H  state). 

It  is  well  known  that  the  Ge  lattice  is  more  flexible  than  the  Si 
lattice,  as  demonstrated  by  its  relatively  lower  bulk  modulus 
(BGe  =  57  GPa  vs.  BSi  =  91.3  GPa)  [31,39].  In  addition,  our  calculation 
shows  that  the  restoring  force  acting  on  a  Ge  atom  displaced  from 
its  equilibrium  position  is  much  smaller  as  compared  to  the  Si  case. 
For  instance,  as  summarized  in  Table  S2,  the  calculated  restoring 
force  for  a  selected  Ge  (Si)  atom  upon  0.02  A  displacements  in  ±x, 
±y  and  ±z  directions  is  around  -0.10  (-0.13)  eV  A'1  in  each  di¬ 
rection.  The  results  provide  clear  evidence  that  the  Ge  lattice  is 
more  flexible  than  the  Si  lattice  thus  able  to  expand  with  less 
associated  resistance. 

Besides  host  lattice  flexibility,  we  find  that  the  Li  bonding 
interaction  varies  substantially  as  the  host  material  changes,  based 
on  the  following  analysis  of  density  of  states  (DOS)  and  charge 
density  distributions.  The  DOS  plot  for  the  Li/Si  system  [Fig.  2(a)  left 
panel]  shows  a  shift  of  the  Fermi  level  above  the  conduction  band 
minimum  of  Si,  indicating  the  charge  transfer  from  Li  to  the  c-Si 
host  matrix.  In  addition,  we  can  notice  significant  overlap  between 
the  Li  2 sp3  and  Si  3sp3  states  in  the  energy  range  of  -5.5  eV 
to  -1.5  eV,  indicating  a  markable  degree  of  covalency  of  the  Li-Si 
bonds.  As  shown  in  Fig.  2(a)  (right  panel),  this  is  also  well  supported 
by  the  isosurface  plot  of  charge  density  differences  (Ap)  before  and 
after  the  Li  insertion;  that  is,  the  charge  accumulation  in  the  region 
between  Li  and  each  of  its  four  nearest  Si  neighbors  with  a  slight 
shift  towards  Si  suggesting  that  the  Li-Si  bonding  has  polar  cova¬ 
lent  character. 

The  DOS  and  Ap  analyses  for  the  Li/Ge  system  are  overall  similar 
to  the  Li/Si  case  but  with  a  few  exceptions  to  distinguish  their 
bonding  properties.  Firstly,  as  shown  in  Fig.  2(b)  left  panel,  the  Ge 
4sp3  peak  is  relatively  narrower  as  compared  to  the  Si  3sp3  case, 
which  is  apparently  due  to  the  relatively  weaker  Ge-Ge  bonding 
interaction  than  the  Si— Si  interaction  as  consistent  with  the  result 
from  the  aforementioned  restoring  force  analysis.  Secondly,  the  Li 
2 sp3  peak  is  noticeably  smaller  in  the  Li/Ge  system  (as  compared  to 
that  in  the  Li/Si  system),  suggesting  the  weaker  hybridization  be¬ 
tween  the  Li  2 sp3  and  Ge  4sp3  orbitals;  this  is  also  consistent  with 
the  smaller  isosurface  volume  of  Ap  with  more  prominent  shift 
towards  Ge  [Fig.  2(b)  right  panel].  The  smaller  degree  of  covalency 
in  the  Li— Ge  bond  relative  to  the  Li-Si  bond  may  allow  Li  to  move 
more  easily;  indeed,  the  restoring  force  experienced  by  a  Li  atom 
displaced  from  its  equilibrium  T-site  by  0.02  A  in  ±x,  ±y  and  ±z 
directions  is  predicted  to  be  about  50%  smaller  than  that  in  c-Si 
[Table  S3]. 

We  further  examined  the  interplay  between  lattice  flexibility  and 
bonding  interaction  by  applying  hydrostatic  strains.  The  lattice 
constant  of  c-Si  was  made  identical  to  that  of  c-Ge,  and  vice  versa, 
causing  5.86%  tensile  and  5.86%  compressive  strains,  respectively. 
Here,  we  will  only  elaborate  on  the  changes  in  tensile-strained  c-Si 
since  the  same  explanation  can  be  applied  to  the  Li/Ge  case  under 
compression.  For  c-Si  under  tension  ( a  =  5.777  A),  as  the  overlap  of 
the  sp3  hybrid  orbitals  decreases  with  elongated  Si-Si  bonds 
[Fig.  2(c)],  so  goes  the  bond  rigidity  as  reflected  by  the  smaller 
restoring  force  (-0.09  eV  A'1  in  x,  y  and  z-directions)  [Table  S2]. 
Likewise,  hybridization  between  the  Li  2 sp3  and  Si  3sp3  states  also 
tends  to  become  weaker  with  the  increased  Li-Si  distance,  as  evi¬ 
denced  by  the  reduced  and  narrowing  Li  2 sp3  DOS  peak  as  well  as  the 
lessening  Ap  isosurface  volume.  As  the  result,  the  restoring  force 
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Fig.  2.  [Left  panel]  Electron  density  of  states  (DOS)  projected  on  Li  and  its  four  nearest  Si  (or  Ge)  neighbors  in  c-Si  (or  c-Ge);  the  intensity  of  Li  2 sp3  DOS  is  scaled  by  4  times.  The 
vertical  dotted  line  indicates  the  Fermi  level  position.  [Right  panel]  Charge  density  difference  (A p)  plots  before  and  after  Li  insertion;  the  red  and  cyan  isosurfaces  represent  the 
regions  of  charge  gain  (+0.0025  e  A-3)  and  loss  (-0.0018  e  A-3),  respectively.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web 
version  of  this  article.) 


experienced  by  the  displaced  Li  atom  (0.02  A  in  ±x,  ±y  and  ±z  di¬ 
rections)  is  found  to  decrease  by  approximately  50%  as  compared  to 
the  strain  free  case  [Table  S3].  While  c-Si  may  become  more  Ge-like 
under  tension,  the  Ea  for  Li  diffusion  is  reduced  to  a  comparable  value 
(0.37  eV)  to  that  in  c-Ge  (0.44  eV).  Contrarily,  the  opposite  effects  are 
true  for  c-Ge  under  compression  (a  =  5.457  A)  [Fig.  2(d)] ;  that  is  c-Ge 
becomes  more  Si-like  due  to  the  increased  Ge-Ge  bond  rigidity  and 
enhanced  Li— Ge  interaction.  The  Ea  value  in  c-Ge  is  found  to  increase 
from  0.44  eV  at  the  strain-free  state  to  0.61  eV  under  compression, 
which  is  remarkably  close  to  the  value  in  c-Si  (0.62  eV).  These 
findings  demonstrate  that  Li  diffusion  is  subject  to  lattice  rigidity  as 
well  as  its  interaction  with  host  atoms;  between  c-Si  and  c-Ge,  a  Li 
atom  can  migrate  more  easily  (lower  Ea)  in  the  later  apparently  due 
to  the  more  flexible  lattice  and  weaker  Li-Ge  interaction. 

3.2.  Dynamic  behavior  in  lithiated  a-Lio.uSi  and  a-Lio.uGe  at  finite 
temperatures 

As  presented  earlier,  our  AIMD  simulations  show  Dy  values  are 
significantly  different  in  Ge  and  Si  during  early  stages  of  lithiation, 


but  with  increasing  Li  contents,  Dy  become  progressively  more 
alike.  For  instance,  Dy  in  a-Lio.i4Ge  is  predicted  to  be 
8.23  x  lO"7  cm2  s-1  at  room  temperature  while  that  in  a-Li0.i4Si  is 
five  orders  of  magnitude  smaller  around  1.60  x  10-12  cm2  s-1.  In 
addition,  we  can  see  that  Dce  is  about  ten  orders  of  magnitude 
larger  than  Dsi,  marking  the  distinctly  different  behaviors  of  a-Si 
and  a-Ge  matrices  particularly  during  early  stages  of  lithiation. 
Therefore,  to  better  understand  the  underlying  reason  for  the 
significantly  enhanced  Dy  and  Dce,  we  analyzed  the  dynamic 
behavior  of  host  lattice  atoms. 

Fig.  3(a)  shows  AIMD  snapshots  of  the  a-Lio.i4Si  and  a-Lio.i4Ge 
systems  annealed  at  800  K;  the  annealing  temperature  was  chosen 
so  that  the  provided  thermal  energy  can  be  sufficient  to  agitate 
atomic  movements  but  not  enough  to  melt  the  lattices.  In  each 
alloy,  we  randomly  selected  a  host  atom  (labeled  as  ‘SU’  and  ‘GeA’  in 
red  color)  and  tracked  its  bonded  neighbors  (marked  as  ‘Sie’  and 
‘Gee’  in  blue  color).  For  a  time  interval  of  8  ps,  in  a-Lio.i4Si,  SU  tends 
to  remain  bonded  to  the  same  four  Sie  neighbors,  suggesting  that 
the  a-Si  host  network  is  nearly  stationary,  barely  rearranges  its 
configuration  in  response  to  moving  Li  atoms.  Contrarily,  GeA  in  a- 
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Fig.  3.  (a)  AIMD  snapshots  of  a-Lio.wSi  and  a-Li0.i4Ge  annealed  at  800  K;  the  host  atom 
of  interest  is  labeled  as  ‘SiA’  and  ‘GeA’  in  red  color,  and  the  host  atoms  involved  in 
atomic  rearrangement  during  a  time  interval  of  8  ps  are  marked  as  ‘SiB’  and  ‘GeB’  in 
blue  color,  (b)  The  extent  of  atomic  rearrangement  (quantified  by  /B)  as  a  function  of 
MD  duration;  here,/B  is  the  averaged  number  of  SiB  (GeB)  neighbors  normalized  with 
respect  to  total  number  of  host  atoms  within  the  supercell.  (For  interpretation  of  the 
references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this 
article.) 


Lio.MGe  is  found  to  undergo  migration  through  a  series  of  bond 
breaking,  bond  switching,  and  new  bond  forming  events.  The 
atomic  rearrangements  associated  with  GeA  extends  far  beyond  the 
initial  adjacent  neighbors,  indicating  a  significant  weakening  of  the 
Ge  lattice  even  at  the  small  degree  of  Li  incorporation,  which  in  turn 
allows  facile  self-diffusion.  To  quantify  the  extent  of  the  atomic 
rearrangements,  the  averaged  number  of  Sie  (Gee)  neighbors 
normalized  with  respect  to  total  number  of  host  atoms  (/b)  within 
the  supercell  is  plotted  as  a  function  of  MD  duration  in  Fig.  3(b); 
these  analyses  were  repeated  multiple  times  on  three  indepen¬ 
dently  constructed  128-atom  supercells.  For  the  system  size  and 
MD  duration  considered  here,  on  average,  a  Ge  atom  in  a-Lio.i4Ge  is 
estimated  to  break  and  form  Ge-Ge  bonds  with  approximately  20% 
of  the  host  atoms  while  that  pertain  to  the  a-Si  host  is  only  around 
4%. 

As  a  result  of  the  facile  atomic  rearrangements,  DGe  is  orders  of 
magnitude  larger  than  DS[  and  consequently  leads  to  faster  Li 
diffusion.  Furthermore,  the  dissimilar  dynamic  behaviors  discussed 
here  may  help  explaining  the  different  mechanical  responses  in  Si 
and  Ge  nanowire  (NW)  lithiation/delithiation  experiments  [24]; 


that  is  (i)  upon  lithiation,  Si  shows  high  anisotropy  with  the  fastest 
and  most  favorable  lithiation  on  {110}  planes  while  Ge  undergoes 
nearly  isotropic  lithiation,  and  (ii)  upon  delithiation,  there  is 
apparent  crack  formation  near  the  center  (along  axial  direction)  of 
the  SiNW  while  the  GeNW  exhibits  porous  structure  instead  of 
cracking.  According  to  our  results,  with  a  small  degree  of  Li  alloying, 
the  rigid  Si  lattice  may  respond  with  less  significant  weakening  as 
compared  to  the  Ge  lattice,  thereby  well  retains  the  crystallo¬ 
graphic  properties  of  unlithiated  planes,  leading  to  a  strong 
orientation  dependence  of  lithiation.  Contrarily,  Ge  shows  pro¬ 
nounced  lattice  weakening  even  at  low  Li  concentrations,  such  that 
upon  lithiation,  the  original  crystallographic  characteristic  is  over¬ 
shadowed  by  the  sufficiently  fast  host  atom  rearrangements, 
thereby  resulting  in  the  isotropic  lithiation.  Likewise,  since  Ge  is 
able  to  easily  undergo  atomic  rearrangements,  the  matrix  is  more 
flexibly  adjusted  to  the  large  strain  variation  during  Li  insertion / 
extraction  and  could  subsequently  reduce  crack  formation. 

Next,  we  looked  into  potentially  important  factors  that  can  in¬ 
fluence  the  dynamic  behavior  of  Ge  and  thereby  Li  atoms.  As 
pointed  out  in  our  previous  study  on  a-LixGe  alloys  [31  ],  the  charge 
transferred  from  Li  can  fill  the  Ge  antibonding  states,  leading  to 
weakened  and  elongated  Ge-Ge  bonds,  which  are  anticipated  to 
have  non-trivial  impacts  on  Ge  atomic  rearrangements;  in  addition, 
the  presence  of  fast  diffusing  Li  atoms  can  also  affect  the  Ge  dy¬ 
namic  behavior.  To  assess  the  relative  contribution  from  each  factor, 
we  calculated  and  compared  DGe  values  under  four  different  sce¬ 
narios  as  listed  in  Fig.  4:  (a)  128-atom  pristine  a- Ge  lattice  with  a 
calculated  mass  density  of  pGe  =  4.87  g  cm-3  for  reference,  (b) 
porous  112-atom  a-Ge  lattice  (pGe  =  4.79  g  cm-3),  simulating  the 
lattice  expansion/weakening  effect  associated  with  elongated  Ge- 
Ge  bonds  (note  that  16  Li  atoms,  which  render  a  lithiated  compo¬ 
sition  of  a-Lio.i4Ge,  are  removed  from  the  system),  (c)  negatively 
charged  112-atom  porous  a-Ge  lattice  (with  16  additional  elec¬ 
trons;  pGe  =  4.79  g  cm-3),  simulating  both  the  lattice  expansion  and 
electron  injection  effects,  and  (d)  128-atom  a-Li0.i4Ge 
(PGe  =  4.79  g  cm-3),  representing  the  full  effects  of  Li  incorporation; 
that  is,  the  Ge  dynamic  behavior  is  also  influenced  by  the  presence 
of  fast  diffusing  Li  atoms  as  well  as  the  contributions  from  lattice 
expansion  and  electron  injection. 

AIMD  simulations  were  performed  at  700  K  to  calculate  DGe  in 
these  four  scenarios.  In  comparison  to  pristine  a-Ge 
(DGe  ~  1.6  x  10-7  cm2  s-1),  the  porous  a-Ge  is  predicted  to  have  a 
higher  DGe  ( ~4.6  x  10-6  cm2  s-1),  indicating  the  lattice  expansion/ 
weakening  effect  tends  to  assist  Ge  self-diffusion  since  the  elon¬ 
gated  Ge-Ge  bonds  are  weaker  thus  able  to  rearrange  more  easily. 
With  the  electron  injection,  the  Ge-Ge  bonds  are  further  weak¬ 
ened,  and  Dee  is  raised  to  7.7  x  10  6  cm2  s-1.  Nevertheless,  DGe  is 
still  much  greater  (~4.0x  10-5  cm2  s-1)  at  the  presence  of  Li  atoms 
in  a-Lio.i4Ge,  suggesting  that  the  facile  rearrangements  of  host 
atoms  are  likely  synergistically  affected  by  fast  diffusing  Li  atoms. 

Lastly,  we  examined  the  effect  of  Si-Ge  alloying  on  Li  diffusivity. 
The  a-Li0.i4Sii_yGey  alloys  of  varying  Ge  contents  y  (generated 
based  on  AIMD  simulations)  are  homogeneous  with  Li  atoms  being 
well  dispersed  in  the  matrices,  which  is  consistent  with  experi¬ 
mental  observations  that  Si  and  Ge  are  miscible  over  the  entire 
composition  range  and  both  are  able  to  form  favorable  alloys  with 
Li  [46,47].  The  predicted  Du,  Ds i,  and  DGe  at  800  K  as  a  function  ofy 
are  shown  in  Fig.  5;  Du  as  well  as  D si  and  DGe  tend  to  rise  mono- 
tonically  with  increasing  y,  exhibiting  a  nearly  linear  relation.  As 
Si-Si  bonds  are  replaced  with  either  relatively  weaker  Si-Ge  or 
Ge-Ge  bonds,  the  host  matrix  becomes  more  flexible  and  easily 
undergoes  atomic  rearrangements,  thereby  resulting  in  the 
enhanced  Li  diffusivity.  While  the  dynamic  behavior  of  the  host 
lattice  is  modified  through  Ge-alloying,  we  anticipate  that  the  facile 
atomic  rearrangement  not  only  can  improve  Li  diffusivity,  but  also 
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Fig.  4.  Schematic  illustration  of  four  different  simulation  scenarios  to  examine  how  self-diffusivity  (DGe)  is  affected  by  the  (b)  lattice  expansion/weakening,  (c)  electron  injection, 
and  (d)  Li-induced  combined  synergistic  effects,  as  compared  to  the  pure  a-Ge  case  [(a)],  together  with  the  predicted  DGe  values  at  700  I<  in  the  lower  panel.  The  total  number  of 
atoms  and  corresponding  Ge  density  (pGe)  for  each  system  are  also  specified. 


Fig.  5.  Predicted  diffusivities  for  Li,  Si,  and  Ge  atoms  in  selected  a-Li0.i4Sii-yGey  alloys 
from  AIMD  simulations  at  800  K. 

contribute  to  better  strain  accommodation,  which  highlights  an 
important  prospect  of  utilizing  Si-Ge  alloy  as  an  anode  material. 
Due  to  the  intrinsically  lower  gravimetric  lithiation  capacity  (of  Ge), 
for  Si-Ge  alloy  anodes,  the  enhanced  Li  diffusivity  is  often  at  the 
expense  of  reduced  initial  capacity.  However,  our  findings  point  out 
that  while  Si-Ge  alloy  anodes  may  have  a  lower  initial  capacity,  the 
improved  strain  accommodation  would  result  in  better  capacity 
retention  and  thus  the  superior  overall  cycling  performance/ 
Coulomb  efficiency;  such  that  beyond  a  certain  number  of  charge 
cycles,  Si-Ge  alloys  anodes  due  to  better  capacity  retention  can 
exhibit  higher  capacity  than  pure  Si  anodes.  That  is  rapid  charging 
rate  and  high-capacity  retention  can  be  achieved  simultaneously 
via  fine-tuning  of  the  Si-Ge  alloying  condition  (composition  ratio 
and  spatial  distribution). 

4.  Conclusion 

We  examined  and  compared  the  dynamic  lithiation  processes  in 
Si  and  Ge  using  DFT  calculations.  The  variations  of  room- 


temperature  Li  diffusivity  (Du)  with  Li  content  (x)  in  a-LixSi  and 
a-LixGe  were  evaluated  using  DFT-based  MD  simulations.  For  Li 
diffusion  in  the  crystalline  matrix,  Du  is  predicted  around 
xlO-13  cm2  s-1  in  c-Si  and  xlO-11  cm2  s-1  in  c-Ge.  With  increasing 
x,  Du  in  a-LixSi  tends  to  rise  by  orders  of  magnitude  from 
xlO-12  cm2  s-1  (x  =  0.14)  to  xlO  7  cm2  s_1  (x  =  3.57),  whereas  Du 
in  a-LixGe  exhibits  virtually  no  concentration  dependence  and  re¬ 
mains  relatively  constant  around  xlO-7.  The  extremely  facile  Li 
diffusion  even  at  low  Li  concentrations  distinguishes  Ge  from  Si  as 
an  anode  material;  this  difference  in  lithiation  behavior  between  a- 
LixGe  and  a-LixSi  becomes  progressively  smaller  with  increasing  Li 
contents.  To  understand  the  underlying  reasons,  we  investigated 
the  Li-host  matrix  interaction  as  well  as  their  dynamic  behavior 
especially  during  early  stages  of  lithiation. 

Firstly,  we  found  that  in  comparison  to  Si,  the  Ge  lattice  is  (i)  less 
rigid  as  Ge-Ge  bonds  are  relatively  weaker/less  directional  than 
Si-Si  bonds,  and  (ii)  interacting  less  strongly  with  Li  since  the  de¬ 
gree  of  covalency  for  Li-Ge  bonds  is  smaller  than  that  for  Li-Si 
bonds;  both  factors  in  turn  allow  Li  to  migrate  more  easily  (lower 
Ea)  in  c-Ge.  Furthermore,  by  applying  5.86%  tensile  (compressive) 
strain  to  c-Si  (c-Ge)  while  c-Si  becomes  more  Ge-like  and  vice  versa, 
the  Ea  values  for  Li  diffusion  also  become  comparable;  Ea  for  c-Si  (c- 
Ge)  under  zero  and  5.86%  tensile  (compressive)  strain  are  0.62 
(0.44)  eV  and  0.37  (0.61)  eV,  respectively.  These  findings  demon¬ 
strate  that  Li  diffusion  is  strongly  subject  to  lattice  rigidity  as  well  as 
its  interaction  with  host  atoms.  Secondly,  in  addition  to  the  superior 
Du  that  is  five  orders  of  magnitude  larger  in  a-Lio.i4Ge  than  a- 
Lio.MSi,  we  found  Dce  to  be  nearly  ten  orders  greater  than  Dsu 
indicating  that  during  early  stages  of  lithiation,  the  a-Si  matrix 
tends  to  be  stationary  relative  to  the  diffusing  Li  atoms  while  the  a- 
Ge  matrix  exhibits  comparably  facile  self-diffusivity.  To  explain  the 
distinctly  different  dynamic  behaviors  of  two  matrices,  a-Lio.i4Ge 
and  a-Li0.i4Si  were  annealed  at  a  slightly  elevated  temperature 
(800  I<)  for  8  ps  to  facilitate  sufficient  atomic  movements.  In 
comparison  to  the  rather  motionless  a-Si  matrix  (each  Si  atom  is 
bonded  with  the  same  four  nearest  neighbors  during  the  entire  MD 
duration),  the  atomic  rearrangements  of  a-Ge  can  extend  far 
beyond  the  initial  adjacent  neighbors  through  a  series  of  bond 
breaking,  bond  switching,  and  new  bond  forming  events,  which  in 
turn  leads  to  facile  self-diffusion  and  consequently  the  larger  Du. 
Later,  we  also  demonstrated  that  the  prompt  rearrangements  of  Ge 
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host  atoms  is  a  synergistic  result  from  the  combined  contributions 
of  (i)  the  weakened  Ge-Ge  bonds/softened  lattice  due  to  the  Li 
incorporation,  and  (ii)  the  agitation  from  the  fast  diffusing  Li  atoms. 
Furthermore,  the  more  flexible  and  facile  rearrangements  of  Ge 
lattice  (as  compared  to  Si  lattice)  also  render  a  probable  explana¬ 
tion  to  the  different  mechanical  responses  in  Si  and  Ge  nanowires; 
for  instance,  anisotropic  SiNW  lithiation  versus  isotropic  GeNW 
lithiation,  and  the  superior  capacity  retention  of  GeNWs. 

Lastly,  we  examined  the  dynamic  behavior  of  a-Li0.i4Sii_yGey 
alloys  to  assess  the  effect  of  Si— Ge  alloying  on  Dy.  While  the  host 
lattice  is  modified  through  Ge-alloying,  Dy  as  well  as  DS[  and  DGe 
are  significantly  enhanced  and  scale  linearly  with  the  increasing  Ge 
content  (y).  This  trend  indicates  that  the  rate  capability  may  be 
appreciably  improved  via  effective  Ge  alloying.  Moreover,  since  the 
Si-Ge  lattice  is  more  flexible  and  can  undergo  faster  rearrange¬ 
ments  upon  lithiation,  we  expect  Si-Ge  anodes  to  outperform  pure 
Si  anodes  in  turns  of  strain  accommodation  upon  lithiation  and 
thereby  the  superior  capacity  retention.  Therefore,  rapid  charging 
rate  and  high-capacity  retention  can  be  achieved  simultaneously 
via  fine-tuning  of  the  alloying  conditions  (composition  ratio  and 
spatial  distribution),  such  that  more  expensive  Ge  is  effectively 
utilized  while  Si-alloying  mitigates  the  cost.  The  origin  and 
extended  impacts  of  the  enhanced  diffusivities  in  Ge  and  Si-Ge 
alloys  are  first  brought  to  light  in  the  present  work,  and  we  antic¬ 
ipate  the  improved  understanding  on  the  lithiation  dynamics  may 
contribute  to  the  design/development  of  the  next-generation,  high- 
power,  and  high-energy  anodes  for  LIBs. 
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